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SCALING  Z-PINCH  PLASMAS  WITH  ATOMIC  NUMBER  FOR  A 
GIVEN  K-SHELL  EMISSION 

I .  INTRODUCTION 

A  major  goal  in  utilizing  z-pinch  plasmas  as  high  brightness 

laboratory  sources  of  x-rays  for  such  diverse  applications  as  x-ray 

lithography,  x-ray  lasers,  materials  studies,  and  plasma  spectroscopy  is  to 

increase  the  radiation  emission  above  about  1  k.eV  through  the  use  of 

moderate  atomic  number  (Z  =  10  to  36)  materials.  Two  kinds  of 

1 —8 

configuration,  single  wires  and  wire  arrays  or  foils  ,  have  been 

investigated  in  the  past.  Gas  puffs  can  also  be  used  in  either  geometry 

though  only  crudely  since  the  supersonic  gas  flow  at  best  approximates 

these  shapes.  In  both  cases,  it  has  been  observed  that  optimal  x-ray 

emission  decreases  on  a  given  machine  as  »^he  atomic  number  of  the  z-pinch 

2  9 

material  is  increased  ’  . 

Single  wire  dynamics  consists  of  a  complicated  interplay  between 
electrical,  MHD,  icnizatioo,  aiivi  radiation  phenomena.  The  question, 
therefore,  of  how  the  x-ray  emission  of  these  materials  will  scale  with 
atomic  number  is,  in  principle,  a  difficult  one.  Wire  arrays,  on  the  other 
hand,  initially  have  a  comparatively  simple  dynamics  (to  zeroth 
approximation).  A  large  fraction  of  the  electrical  energy  is  first 
converted  into  the  kinetic  energy  of  implosion,  which  is  then  thermalized 
when  the  array  collides  with  itself  at  the  axis  of  the  pinch.  In  this 
case,  a  force  and  a  circuit  equation  suffice  to  describe  the  optimal 
electrical  coupling  of  the  z-pinch  load  to  the  pulse  power  machine^ In 
this  paper,  we  will  describe  how  the  model  of  Ref.  10  can  be  modified  and 
extended  to  provide  predictions  of  how  machine  and  diode  designs  must  be 
scaled  in  order  to  maintain  a  given  K-shell  x-ray  emission  as  the  plasma 
atomic  number  increases  in  z-pinch  array  implosions. 

II.  THE  MODEL 

Let  m  denote  the  ma'-s  -'er  unit  length  of  the  wire  array  or  foil  and 
let  r(t)  be  the  outside  rc..iUS  of  the  circular  array  as  a  function  of  time. 
In  response  to  the  magnetic  forces  that  are  generated  when  a  current  I(t) 
flows  through  the  array,  it  accelerates  inward: 

Manuscript  approved  April  26,  1989. 
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d^r  1  3l'  2 

^2  =  2  3r 


(1) 
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where  L'  is  the  array  load  inductance  per  unit  length  .  For  an  array 


consisting  of  N  wires 


Ik 

3r 


N-1  %  1 


N  2it  r 

Following  Ref.  10,  we  will  utilize  the  dimensionless  variables, 


X  =  r/r  , 
o 


-t  =  t/t  , 
o 


i  =  I/I  , 
o 


where  r^  is  the  initial  array  radius  and 


N-1  o  0 


N  2nR  R 
o  0 


(2) 


(3) 


(4) 


is  the  maximum  current  that  flows  in  the  load,  t  is  the  length  of  the 
array  load,  R^  is  the  diode  resistance,  and  s  p^£/(2n).  In  terms  of 
these  variables,  Eq.  (1)  becomes 


dix 

dx^ 


.2 

-  A  i- 
x 


(5) 


where,  for  large  N,  the  dimensionless  constant  A  is  given  by 


A  = 


2<V^>^o 

m(r  /t 
'  o  o' 


(6) 


For  a  vide  variety  of  existing  pulse  power  machines,  the  current  pulse 

0 

is  adequately  described  by  a  linear  ramp  I(t)  =  I^t  during  the  main  run-in 
phase  of  the  array  implosion.  In  pulse  power  machines  that  are  designed  as 
current  sources,  this  should  also  be  a  useful  description  of  the  current 
pulse.  From  Eq.  (6),  one  can  derive  an  equation  for  the  increase  of 
kinetic  energy  with  implosion  distance: 


2 


1 

2 


(7) 


It  is  clear  from  Eq.  (7)  that  maximum  kinetic  energy  is  generated  in  the 
implosion  ir  the  current  is  large  in  the  final  stages  of  the  run-in.  Thus, 
on  both  experimental  and  theoretical  grounds,  it  is  useful  to  investigate 
the  scaling  of  z  pinch  implosions  for  linearly  rising  current  pulses.  For 
this  case,  one  can  define  a  new  dimensionless  quantity  q  by 

q  =  Qq  t  ,  (8) 

where 


.1/4  °  ^  ,1/2 

q  s  A  (I  t  /I  ) 

O  0  0  0 


(9) 


so  that  Eq.  (5)  simplifies  to 
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d^x 

dq^ 


(10) 


The  picture  of  array  implosions  that  is  provided  by  Eq.  (10) 
eventually  breaks  down  when  the  plasma  begins  to  collide  on  itself  and  to 
thermalize  the  kinetic  energy  generated  prior  to  this  collision. 
Experimental  evidence^*^  on  implosion  times  obtained  from  a  number  of 
machines  suggests,  in  fact,  that  Eq.  (10)  is  generally  valid  until  the 
plasma  has  collapsed  to  about  10%  of  its  initial  radius,  at  which  time 
further  collapse  is  no  longer  discernible  and  the  main  radiation  pulse 
commences.  Thus,  we  will  be  interested  in  the  solution  to  Eq.  (10)  up  to 
values  of  x  -  0.1  subject  to  the  boundary  conditions  that  the  implosion 
begins  at  rest  dx/dT  =  0  when  x  =  1. 

Eq.  (10)  can  be  solved  either  analytically  by  power  series  or 
numerically.  The  solution  information  of  interest  is  the  final  values  of  q 
and  (dx/dq),  q^  and  (dx/dq)^  respectively,  that  are  achieved  when  x  reaches 
its  final  value  x^  just  prior  to  the  thermalization  process.  These 
quantities,  obtained  from  numerical  solutions  to  Eq.  (10),  are  shown  in 
Figures  1  and  2.  Least  squares  fits  to  the  computed  values  of  q^  and 
(dx/dq)j  are  also  drawn  in  these  figures.  These  fits  are  given  by  the 


1 


3 


polynomial  expressions, 


=  1.73  -  0.22  x^-  0.42  x^^  ,  (11) 

(dx/dq)j  =  4.72  -  18.9  Xj+  49.1  ,  •  (12) 

valid  for  0.04  <  x^  <  0.14. 

III.  Scaling  Arguments 

By  assuming  a  linear  rise  in  current,  one  guarantees  that  the  time 
of  array  implosion  t^  will  coincide  with  the  time  of  peak,  current.  What  we 
would  like  to  require,  in  addition,  is  that  during  the  implosion,  the  array 
acquires  enough  kinetic  energy  to  ignite  the  K-shell  and  contains  enough 
mass  to  produce  the  K-shell  emission  rates  that  are  required  of  a  high 
yield  super-kilovolt  x-ray  source.  In  the  process,  we  would  like  to 
determine  how  these  requirements  scale  with  atomic  number,  Z,  from  Z  =  10 
(neon)  to  Z  =  36  (krypton). 

To  begin,  consider  a  uniform  plasma  cylinder  of  unit  length  and 
radius  r^  containing  a  single  element  of  atomic  number  Z.  The  amount  of 
mass  m  contained  within  this  cylinder  is 

ra  =  Am^  Nj  It  rp  (13) 

where  A  is  the  atomic  weight  of  the  element,  m^  is  the  mass  of  a  proton, 
and  is  the  ion  density  of  the  plasma.  A  least  squares  polynomial  fit  to 
log  A  as  a  function  of  Z  results  in 

A  =  1.58  Z  (14) 

This  functional  fit,  up  to  Z  =  36,  to  the  atomic  weights  is  shown  in  Figure 
3 .  Thus 

m(pg/cm)  =  8.28x10  r  ^Z  (15) 

In  this  paper,  we  will  employ  a  mixed  set  of  MKS  units  that  are  appropriate 
to  a  discussion  of  2-pinches.  Thus,  we  will  generally  express  distances 
and  ion  densities  in  units  of  centimeters  and  cm  ,  masses  in  units  of 


micrograms,  yields  in  kilojoules,  energies  and  temperatures  in  eV,  currents 

in  mega-amperes,  and  time  in  nanoseconds  or  microseconds. 

There  are  two  ways  to  determine  the  scaling  of  K-shell  yield  with 

atomic  number.  One  is  by  a  simple  theoretical  argument,  the  other  is  by  an 

analysis  of  previous  and  future  experimental  and  theoretical  data.-  Because 

2 

of  its  large  oscillator  strength,  the  He-a  Is  -ls2p  line  dominates  the  K- 

shell  emission  from  low  density,  optically  thin  plasmas.  Thus,  one  can 

argue  on  the  basis  of  the  optically  thin  emission  rate  of  this  line  that 

the  yield  per  unit  length,  y,  of  total  K-shell  emission  will  scale  as  a 

power  of  and  inversely  as  a  power  of  Z.  The  coefficient  of 

proportionality  can  be  determined  either  from  direct  spectroscopic 

measurements  of  N.  and  T  made  in  previous  low  atomic  number  z-pinch  shots 
1  e  . 

for  which  an  optimal  K-shell  yield  was  approximately  obtained  or  from 
theoretical  calculations.  One  then  finds  that 


y(kJ/cm)  =  Cy  n“  Z  ^(Z)  rj 


(16) 


where  E(Z)  is  an  exponential  function  of  Z  to  be  defined  shortly. 

Estimates  for  the  exponents,  a,  0,  and  y  can  be  obtained  from  the 
following  coronal  equilibrium  arguments.  As  the  atomic  number  of  a  given 
sized  plasma  is  increased,  the  photon  probabilities  of  escape  in  the 
resonance  lines  of  helium  or  hydrogen  approach  1.  Thus,  the  K-shell 
radiation  rates  of  the  plasma  become  increasingly  optically  thin.  In  this 
case,  the  amount  of  He-a  emission  from  our  uniform  plasma  cylinder  of  unit 
length  will  be 


^He-a  =  ^He  ^  ^ls-2p  ^  ""f  \ 


(17) 


where  is  the  number  of  ions  in  the  helium-liko  ground  state,  is  the 
electron  density,  X2g_2p  rate  coefficient  for  exciting  the  ls2p^P 
state  from  the  ground  state,  AE  is  the  energy  of  the  transition  and  of  the 
radiated  photon,  and  is  the  duration  of  the  emission.  To  determine  the 
scaling  of  y„^  with  N.  and  Z  we  use  Seaton's  formula  for  the  collisional 
excitation  rate  for  dipole  allowed  transitions^  : 


-AE/T 

X.  -  (cm'^/sec)  =  1.6x10  fg  ^ ,  (18) 

is-zp  ^ 

e 

where  f  is  the  oscillator  strength,  g  is  the  thermally  averaged  gaunt 
factor  ~  1,  and  T  is  in  units  of  eV.  The  temperature  T  for  which  K-shell 

S  0  o 

occupation  and  excitation  is  strong  does  not  scale,  as  AE  does,  with  Z" 
because  of  the  different  Z  dependences  of  the  collisional  ionization  and 
radiation  recombination  rates.  Based  on  estimates  of  ionization 

equilibrium  in  hydrogen,  we  will  take  this  scaling  to  be 


Tg(eV)  =  0.3  2  (19) 

although  a  more  detailed  analysis  of  this  scaling  should  be  made.  Thus, 

Xis_2p  '  E(Z)  /  (20) 

where,  because  AE  ~  Z^  and  T^  -  Z^'^,  we  take  E(2)  =  exp(-20.6/Z*^’'^) .  In 

order  to  determine  the  scaling  of  with  Z,  we  assume  that  K-shell 

emission  is  proportional  to  the  inertial  confinement  time  of  the  plasma 
during  emission; 


f 


K  v^,  /-TIT.  /.  vl/2 

th  (2kT  /Am  ) 

e  p' 


-  r^/  Z 


0.9 


(21) 


Thus,  in  coronal  approximation. 


Nfr^  E(Z) 
^He-a  ~  “^1^35 


(22) 


and,  in  this  case,  taking  y  and  scale  identically  with  Z  and 

we  find  that  a  =  2,  6  =  1.35,  and  y  =  3. 


If  the  K-shell  of  a  given  element  is  to  be  excited,  a  minim.um  energy 
must  be  imparted  to  each  imploding  ion  in  order  to  (1)  strip  all  but  one  or 
two  electrons  from  the  ion  and  (2)  maintain  each  stripped  electron  at 
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temperatures  to  avoid  recombination  back  into  the  L-shell.  Of  course 
i nef f i r i enc i es  in  the  r herma H za r i on  process,  in  particular  subkilovolt 
radiative  losses,  imply  that  much  more  energy  than  this  is  needed  in 
practice.  However,  an  initial  insight  into  basic  Z-scalir.g  may  be  obtained 
by  consideration  of  these  basic  quantities. 

As  noted,  the  temperatures  at  which  maximum  amounts  of  K-shell 

9  Q 

radiation  are  expected  are  =0.3  2“" ' ■  We  approximate  the  number  ot 

stripped  electrons  as  equalling  2.  Therefore,  the  thermal  energy  of  the 
stripped  electrons  from  an  ion  of  atomic  number  Z  at  this  temperature  is, 

Ej.  =  0,4  Z^'^  eV/ion.  (23) 

The  energy  needed  to  strip  the  electrons  from  each  ion  is  also  important 
and  non-negligible.  We  shall  assume  that  the  plasma  is  50%  hydrogen- like 

and  50%  helium-like.  The  ionization  potentials  for  each  element  are 

1 1, 

available  from  many  references;  we  have  used  the  data  of  Carlson  et  al. 
The  energy  needed  is  given  by  the  following  fit  -  which  for  purposes  of  the 
approximate  treatment  of  this  paper  has  negligible  ertv^r; 

E  =  3.83  Z^'^^  eV/ion.  (24) 

s 

Obviously,  the  minimum  energy  per  ion  that  must  be  acquired  during 
thermalization  is  given  by  E^+E^.  A  critical  question  is,  in  practice 
however,  how  much  energy  per  ion  beyond  the  minimum  is  required  for  optimum 
K-shell  radiation  to  be  generated  in  a  z-pinch?  We  will  assume  that  the 
kinetic  energy  per  ion  that  must  be  generated  in  an  array  implosion  is  some 
multiple,  E,  of  this  minimum  energy.  From  past  analyses  of  aluminum 
experimental  data,  there  is  some  indication  that  e  -  15  and  that  it  may  be 
as  large  as  30.  We  now  make  the  assumption  that,  for  all  Z,  e  (~  15)  times 
(E^+E,)  will  be  required  for  efficient  K-shell  radiation  production.  Such 
a  large  factor  is  plausible  when  one  considers  the  numerous  sinks  for 
generator  energy  other  than  production  of  K-shell  radiation.  Such 
processes  include  runaway  electron  producrinp  during  the  implnp’'»^  fhe 
generation  of  axial  plasma  motions,  nonuniform  thermalization  and 
radiative  losses  other  than  K-shell.  Thi<;  is  not  unreasonable  as  a  first 
approximation  since  it  allows  radiative  and  other  losses  to  scale  upward 
with  Z  in  the  same  fashion  as  the  thermal  and  ionization  energies,  which  is 
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initially  rpa.^o'ipbl  e .  Beginning  vith  a  factor  of  15  can  be  viewed 
optimistically,  since  it  indicates  that  improvements  in  the  efficiency  of 
these  devices  may  be  achieveable.  This  assumption  means  that 

1  2  -19  19  7  hU 

^  m.  (dr  .t)p  =  (KE).  =  1.6  x  10  (0.4  3.83  ®  ) 

2  1  t  ion 

=  2.39  X  10"^'  E  Z^’^^  (J/ion),  (25) 

where,  in  Eq.  (25),  the  additional  simplification  of  the  power-law  sum  to  a 
single  power  law  does  not  result  in  substantial  error  for  present  purposes. 

IV.  RESULTS 

Equations  (11),  (12),  (14),  (15),  (16),  (19),  and  (25)  allow  one  to 

determine  the  scaling  of  the  mass  loading,  implosion  time,  and  the  rate  of 

O 

current  rise,  (I-dot),  with  atomic  number  in  order  to  achieve  a  given  K- 
shell  yield.  For  example,  on  eliminating  from  Eqs.  (15)  and  (16),  one 
finds  that 


m  (ng/cm) 


=  C 


1/a 


(2-3/a)„ 
f  ^ 


(l.U0/a)g(Z) 


-l/o 

» 


(26) 


where  C  =  8.28xlO~^^/C  The  implosion  time  is  t,  =  (t  /q  )  9,;  it  can 

m  y  f  o  0  f 

be  found  from  Eqs  (11),  (12),  and  (25)  since 


Thus , 


(dr/dt)  j 


(dx/dq)^ 

o 


— —  (dx/dq),. 


(27) 


f  —  1  ( 

,  m. 

1  dq  1 

2(KE).^J  -^0 

where  we  now  write  (dx/dq)^  for  the  magnitude  of  this  final  velocity. 
Finally,  we  note  that  the  I-dot  of  the  implosion  can  be  dete*rmined  from 
Eqs.  (6),  (9),  (28)  and  the  relationship  t./qc  =  t  /q  : 

L  L  O  O 
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For  a  given  yield  and  initial  array  radius,  all  of  the  quantities  on 
the  right  hand  sides  of  Eqs-  (26),  (28),  and  (30)  are  scalable  with  Z, 
assuming  that  conventional  wisdom  is  correct  and  back  EMF's  and  plasma 
pressure  forces  will  stop  the  1-d  implosion  and  begin  to  thermalize  the 
kinetic  energy  at  or  around  a  value  of  -  0.1  (see  Ref.  10).  In  order  to 
illustrate  this  procedure  in  more  detail,  ve  will  now  specialize  to  two 
cases  where  a=2,  3=1.35,  and  y=3  (see  Eq.  22)  and,  from  experimental 
estimates,  we  take  to  be  9.1x10“^^.  Eqs.  (26),  (28),  and  (30)  then 
become 


m  (ug/cm)  =  8.68  x  10  "^(yr)^'"^  exp(9 . 2/Z°’ ^) ,  (31) 

t^(ns)  =  7.5  X  10^  q^  (dx/dq)^  (32) 

,  ,1/4  ,3.3 

o  -5  2  ^  0  Q 

I^(MA/ms)  =  4.44x10  (dx/dq)j  e  - ^ -  exp(4.6/Z  ),  (33) 

0 

where  as  mentioned  earlier,  y  is  expressed  in  units  of  kJ/cm  and  r^  and  r^ 
in  units  of  cm.  Graphs  of  these  scaling  relations  are  shown  in  Figures  4-6 
for  two  cases. 

In  case  I,  we  took  the  initial  radius  of  the  array  to  be  1  cm,  and 
assumed  Xj=0.1,  so  that  qj=1.7,  (dx/dq) ^=3.3,  and  r^=0.1  cm.  The  desired 
yield  was  taken  to  be  2kJ/cm  and  the  thermalizat ion  efficiency,  e=15.  In 
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case  II,  these  values  veie  taken  to  be;  r  =2  cm,  Xc=n.07,  q,=1.71, 

O  L  L 

(vlx/dq )  ^-  =  3 . 63 ,  r^=0.14  cm.  y=15  kJ/cm,  and  e=10. 

It  is  clear  from  Eqs.  (31)  -  (33)  and  <^Lom  Figs.  4-6  that  array 
implosions  must  start  at  a  large  radius  and  implode  to  a  fairly  tight 
cylinder  on  axis  with  a  gooa  kinetic-to-thermal  conversion  efficiency  in 
order  to  produce  high  K-shell  yields  with  long  implosion  times  and  small  I- 
dots.  Moreover,  for  a  given  machine,  with  given  I-dot  and  current 
capabilities,  the  optimum  choice  of  element  depends  on  the  yield-per-length 
and  spectrum  requirements.  It  is  also  clear,  however,  that  the  pulse  power 
machine  must  be  redesigned  as  the  load  element  or  yield  requirement  is 
changed. 

The  above  scaling  relations  derive  from  the  assumptions  of  an 

optically  thin  plasma  and  an  impressed  linearly  rising  current  across  the 

plasma  that  terminates  at  the  time  of  kiretic  energy  thermalization. 

Figure  7,  which  shows  the  ion  densities  that  are  computed  from  Eq.  (15)  and 

(31),  indicates  that  the  first  assumption  holds  up  in  the  two  cases  under 

study  very  well.  As  to  the  second  assumption,  we  replaced  the  voltage 

1C 

source  utilized  in  Katzenstein's  work  with  a  current  source  with  obvious 
fundamental  consequences,  since  array  implosions  are  current  driven. 

In  general,  the  ohmic  r'-sistance  of  the  array  plasma  is  small,  so 
that  the  voltage  drop  across  the  load  V^(t)  is  given  by  =  d(L^I)/dt, 

wh*-  the  load  inductance  L.  for  an  N-wire  array  of  wire  diameter  d  having 

^10 
a  return  current  path  radius  R  is  given  by  : 


Hence 


and 


=  L 

f  J-1 

N-1  , 

+  ”7; —  In 

f  -  11 

0 

(  N 

1  Nd^J 

N 

l  r  )) 

^  -  L  i  f- 

dt  ~  o  r  I  dtj  ’ 

=  I^  (t/r)(-dr/dt)). 


(34) 


(35) 


(36) 


After  converting  to  dimensionless  variables,  one  finds  that  the  voltage 
drop  across  the  load  at  the  end  of  the  implosion,  V„j,  is 
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=  (4  +  (qf/Xf)(dx/dq)£  lJ  .  (37) 

3  3 

By  virtue  of  Eq.  (33),  therefore,  it  follows  that  V  ,  scales  as  Z  .  . 

In  general,  the  second  term  in  Eq.  (37)  dominates  over  the  first.  As 

the  load  collapses,  therefore,  maintenance  of  a  steadily  rising  current 

requires  a  disproportionately  large  rise  in  the  voltage  across  the  plasma 

near  the  end  of  the  collapse.  This  point  is  demonstrated  in  Figures  8  and 

0 

9.  We  have  defined  the  peak  current  to  be  For  an  array  length  of 

2  cm,  we  see  from  Figures  8  and  9  that  peak  currents  of  near  2  MA  in 

aluminum  generate  peak  voltages  at  the  end  of  the  collapse  near  7  MV.  For 

selenium,  however,  peak  currents  of  11  MA  and  14  MA  are  needed  in  the  two 

cases  and  they  produce  peak  voltages  near  150  MV.  Under  these  conditions, 

pulse  power  machines  will  not  act  as  current  sources  and  the  current  will 

generally  take  a  steep  dip  at  the  time  of  the  final  collapse  of  the  plasma 

on  axis,  rendering  the  assumption  of  a  linear  current  rise  invalid.  For 

this  reason,  the  scaling  relations  developed  in  this  paper  will  generally 

need  some  correcting  and/or  to  be  empirically  determined  under  experimental 

conditions  that  implode  the  arrays  in  the  required  times  and  give  them  the 

required  final  kinetic  energy  per  ion. 

One  theoretical  feature  of  the  Eqs.  (19)  -  (33)  scaling  relations, 

already  has  some  experimental  support:  that  of  the  scaling  of  yield  with 

0 

peak  current  for  optimal  yields.  Since  I^=I^tp  it  follow  from  Eq  (30) 
that 


Therefore, 

that 


m(kg/m)  =  ^ 


2  I 


(38) 


from  Eqs.  (13)  and  (16)  (with  (x=2,  0=1.35,  and  y=3)  one  finds 


y(kJ/cm)  = 


2.6x10  E(Z) 


rfZ 


8.37 


(39) 
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Thus,  the  scaling  relations  can  be  rewritten  to  exhibit  an  scaling  of 

the  yield  on  peak  current.  In  the  context  of  a  fixed  pulse  power  machine, 

such  scaling  has  been  observed  experimentally  for  optimal  K-shell  yields  in 
14 

neon. 

It  is  also  of  interest  to  determine  what  the  scalings  of  Eq.  (31)  - 
(33)  predict  on  average  about  the  balance  of  pressure,  f^,  and  magnetic 
field,  forces  per  unit  length  at  the  end  of  the  implosion.  Using 

Eqs.  (13)  and  (19),  we  find  (in  HKS  units)  that 

f  =  2nr,p,  =  2iir,ZN.kT 
p  f'^f  fie 

-2n  Q 

=  4.8x10  (2m/(m.r^))2^‘^.  (40) 

Also,  from  Eqs.  (25),  (28),  and  (30), 


f  -  ^ 

^jxB  “  4rtr^ 

=  2.39x10"^^  f  — 

Ufr^ 

These  forces  will  be  equal  only  for  a  particular  thermalization  efficiency 

2 

e  =  0.2  [dx/dq]  q^,  (42) 

which,  in  the  two  cases  under  consideration  in  this  paper,  has  the  values 
(<4)  displayed  in  Figure  10.  Since  E  is  expected  to  have  values  much 
larger  than  those  in  Figure  10,  f.  will  be  much  larger  than  f  .  However, 

jxd  P 

as  noted  above,  in  reality,  the  current  will  dip  sharply  at  the  moment  the 
plasma  arrives  on  axis.  Thus,  the  collapse  will  be  stopped  as  the  forces 
reverse  in  strength.  This  accounts,  to  some  extent,  for  the  empirical 
observation  that  x^-0.1.  Plasma  instability  development  may  also 
contribute  to  this  effect. 
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V.  Discussion  and  Conclusions 


The  basic  idea  underlying  the  scaling  relations  of  Eqs.  (31)  -  (33) 
is  that  if  z-pinch  x-ray  sources  with  a  given  yield  and  successively  higher 
kilovolt  energies  are  to  be  developed,  then  one  must  redesign  -the  pulse 
power  machine  as  the  atomic  number  of  the  plasma  is  increased  in  order  to 
implode  the  required  higher  plasma  mass  with  a  higher  power  input  and  a 
shorter  implosion  time  (for  a  given  plasma  load  geometry).  It  follows, 
conversely,  that  if  the  machine  is  not  redesigned  or  reconfigured  to  the 
plasma,  that  one  can  expect  the  yield  of  K-shell  x-rays  to  fall  as  the 
atomic  number  of  the  plasma  is  increased  in  conformity  with  experimental 

9 

observations.  It  also  follows  that  the  theoretical  convenience  of  a 
linearly  rising  current  pulse  should  not  be  strictly  required 
experimentally,  provided  the  required  mass  can  be  imploded  to  the  required 
final  kinetic  energy  by  an  appropriate  circuit  design. 

The  more  fundamental  scaling  relations  of  this  paper  are  given  by 
Eqs.  (26),  (28),  and  (30)  with  a,  8,  y,  c,  x^,  q^,  and  (dx/dq)^  regarded  as 
empirical  constants  to  be  determined  ultimately,  either  from  experimental 
data  or  from  1-d  hydrodynamic  computer  simulations  of  the  array  or  foil 
implosion  process.  Simulations  are  needed,  in  particular,  in  order  to 
determine  (1)  the  thermalization  efficiency,  e,  of  foil  implosions  and  (2) 
whether  or  not  this  efficiency  itself  is  a  scalable  function  of  Z  and/or  a 
function  of  yield.  We  carried  out  some  preliminary  1-d  calculations  for  a 
titanium  foil.  They  suggest  that  thermalization  efficiencies  of  15  do 
occur  for  a  certain  class  of  implosions.  One  major  limitation  of  these 
calculations  is  that  they  do  not  describe  the  instability  growth  that  is 
known  to  occur  experimentally,  and  that  undoubtedly  plays  a  role  in 
limiting  the  application  of  the  above  theory  and  the  extent  of  the 
implosion  (i.e.  the  size  of  x^)  and,  thus,  the  thermalization  efficiency. 
However,  it  should  also  be  noted  that  experiments  conducted  to  date,  in 
which  instability  growth  was  observed  and  in  which  the  implosion  was 
limited,  were  generally  non-optimally  designed  (as  described  in  this  paper) 
for  radiation  production  in  the  K-shell.  Another  factor  that  may  limit  the 
range  of  validity  of  these  scaling  relations  is  the  tendency  of  z-pinches 
to  undergo  radiative  collapse. If  the  mass  loading  is  sufficiently 
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high  or  if  instabilities  form,  the  plasma  will  implode  to  radii  much 

smaller  than  D.l  r  . 

o 

Finally,  important  experimental  evidence  of  non-optimum  performance 

is  frequently  provided  by  filtered  pinhole  pictures  of  the  z-pinch  emission 
2  3 

region.  ’  It  is  generally  observed  that  as  the  energy  of  the  photons 
increases,  the  size  of  the  corresponding  emission  region  decreases. 
Depending  on  the  pulse  power  machine,  this  generally  means  that  as  Z 
increases,  the  fraction  of  the  foil  mass  that  participates  in  the  K-shell 
emission  decreases.  In  this  sense,  Z-pinch  plasmas,  like  vacuum  sparks, 
are  remarkably  accommodating.  When  more  plasma  than  can  be  ionized  to  the 
K-shell  is  imploded,  only  an  appropriate  fraction,  matched  to  the  power 
capabilities  of  the  circuit,  is  heated  to  produce  the  K-shell  emission. 
This  experimental  behavior  contrasts  with  one  of  the  basic  assumptions  of 
the  model  calculations  of  this  paper  that  all  of  the  foil  mass  that  is 
accelerated  inward  participates  in  the  subsequent  thermalization  and  K- 
shell  emission.  Any  experiments  that  are  conducted  to  test  the  scaling 
r'^lations  proposed  in  this  paper  should  carefully  examine  the  validity  of 
this  mass  participation  assumption. 
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Figure  2.  A  least  squares  polynomial  fit  to  calculated  values  of  (dx/dq)^ 
is  shown  valid  for  values  of  between  0.04  and  0.14. 
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Figure  3.  The  atomic  weights  between  neon  (Z=10)  and  krypton  (Z=36)  are 
drawn  relative  to  the  curve  fit  A-1.58Z^’^. 


Figure  4.  The  scaling  of  the  mass-per-length  of  a  z-pinch  array  as  a 
function  of  atomic  number  is  shown  for  the  two  cases  (I  and  II)  described 
in  the  text. 
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ligiiie  7.  riie  final  ion  density  tliat 
i  nip  1  os  i  oin: ,  i  and  II,  is  shown  as  a  fiinctio 
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Figure  10.  The  thermalization  efficiency  that  is  needed  to  achieve  force 
balance  at  the  final  radius  of  the  implosion  is  shown  as  function  of  atomic 
number,  for  the  two  cases,  I  and  II. 
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